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= COMPLEX*8 ** COMPLEX*B power routine 


DECLARAT JONS 
OTSSPOWCC = COMPLEX*B8 ** COMPLEX*8 power routine 
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TITLE OTSSPOWCC = COMPLEX*8 *# COMPLER NG ower routine 
ENT /1-004/ 3; File: OTSPOWCC.MAR 


SRRAAAAHSS AHHH TAAAAALE AAA AAAAAAAAAAAAAAAAAAAAAAARAERERAAAAAREAAAAAAAAAAAEee 


COPYRIGHT («) 1978, 1980, 1982, 1984 BY 
UIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 


te ® 

;* a 

z* DIGITA L FO * 

, :* ALL RIGHTS RESERVED. . 

5 2 

11 ;* THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED * 

12 ;* ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE ®* 

15 ;* INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER * 

14 ;* COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY * 

° 15 ;* OTHER PERSON. NO TITLE TO wn OWNERSHIP OF THE SOFTWARE IS WHEREBY * 

‘4 ! :* TRANSFERRED. e 

g * 

00 18 ;* THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 

00 19 ;* D SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT * 

Bh :* CORPORATION. * 

3 o 

0 3* DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS * 

BS ” SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DiGITAL. . 
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++ 
; FACILITY: Language support Library - user callable 


; ABSTRACT: 
COMPLEX*8 base to COMPLEX*8 power giving COMPLEX*8 result. 
; ENVIRONMENT: User Mode, AST Reentrant 
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; MODIFIED BY: 
; i=) pry BH VERSION 00 


So 


. le 

. $2 - Statdordized version number format, with three digits in 
: the edit cy! JBS 16-NOV- 

; 1-003 - add "to the PSECT directive. JBS oe DeC~28 

; 1-004 - Changé shared external references to G* RNH 25-Sep-81 
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: AUTHOR: Steven B. Lionel, CREATION DATE: 24-Oct-1978: Version 0 
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. -SBTTL OTSSPOWCC = COMPLEX*8 ** COMPLEX®*8 power routine 
FUNCTIONAL DESCRIPTION: 

OTSSPOWCC evaluates the result of taking a gone lex base 

to a complex power. The ANS FORTRAN X3.9-1978 standard defines 
complex exponentiation as: 

x #@ y = CEXP( y © CLOG(x)) 

where x and y are type COMPLEX. 


] 
| 
The arguments of OTSSPOWCC are CALL BY VALUE. 
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CALLING SEQUENCE: 
power.wfc.v = OTSSPOWCC (base.rfc.v, exponent.rfc.v) 
INPUT PARAMETERS: 
Both base and exponent are COMPLEX*8 numbers, each consisting 
of a REAL*4 real part and a REAL*4 imaginary part. Both are 
CALL BY VALUE. 
IMPLICIT INPUTS: 
NONE 
OUTPUT PARAMETERS: 
NONE 
IMPLICIT OUTPUTS: 
NONE 
FUNCTION VALUE: 


The COMPLEX*8 (REAL*4, REAL*4) result of taking the 
COMPLEX base to the COMPLEX exponent power. 


SIDE EFFECTS: 
Possible error signals are: 


MTHS_INVARGMAT if base is (0.,0.). 

MTHS_FLOOVEMAT if floating overflow occurs. 

MTHS_SINCOSSIG if absolute value of the imaginary part of 
(exponent * CLOG(base)) > PIl#2**30. 

$S$_ROPRAND if reserved floating operand is fetched. 
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sENTRY OTSSPOWCC, “M<R2,R3,R4,R5> ; disable integer ovflo 
MTHSFLAG_JACKET 3 establish math error handler 
MOVAB G*MTHSSJACKET_HND, (FP) 


set handler address to jacket 
handler 


mMOVa exp(AP), R2 ; put exponent in R2, R 
ater operations will 
for reserved operands 


3 
check 


Get complex logarithm of base 


PUSHAL base(AP) ; address of base 
CALLS #1, G*MTHSCLOG ; RO,R1 get LOG(base) 
; call by reference 


ye? eae is in RC, R1. Multiply by exponent. 
RO,R1 = CLOG(base) = a+bi 

R2,R3 = exp = c+di 

Complex multiplication defined as: 


real part = ac-bd 
imaginary part = ad+bc 


MULF3 RO, R3, R4 ; R4 = ad 
MULF Re, Re ; RO = ac 
MULF R1, R35, RS : RS = bd 
SUBF RS, RO ; RO = ac-bd 
MULF R2, R1 ; RI = be 
ADDF R4, R1 ; R1 = adtbc 


Now compute CEXP(product) 


=(SP) ; put product (RO,R1) on stack 


mOVa RO, 3 

PUSHL SP 3 agers s of ergynents 

CALLS #1, G*MTHSCEXP ; RO, R1 get EXP(product) 
g call by reference 

RET 3; all done, exit 
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Symbol table Ate aes gree tEs 1i3Fi28 HAXANES Bacte VOR OO ans1 POM 8) 1- 


BASE = 

EXP = Cc 
MTIHSSJACKET_HND teeeeere 
MTHSCEXP 


MTHSCLOG teeeeere G 
OTS$SPOWCC 00000000 RG 
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PSECT name Allocation PSECT No. Attributes 


- ABS, 444! ( 0.) 60 ( Q.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
-OTSSCODE 00000 C $6.) OFC 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 


Gewoccccccceccccccccccccns 


Phase Page faults CPU Time Elapsed Time 


Initialization % 80:00:08 .08 00:00:01.80 
Command processing 1 00:00:00.59 00:00:05.14 
Syabol table sort % §0:00:00:60 $0 09;00:00 
r $ $ e : + e 

Ri 00:00:00-48 00:00:04.60 
ed table output 00:00:00.01 00:00:00.05 
Psect synopsis output 00: 80005 00:00:00.03 
Cross~reference output 3 00:00:00.0 00:00:00.00 
Assembler run totals 30 00:00:01.85 00:00:15.34 


The_ working set Limit was 900 pages. 

2597 bytes (6 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space allocased to hold 6 non-local and 0 Local symbols. 
235 source Lines were read in Pass 1, producing 11 object records in Pass 2. 

1 page of virtual memory was used to define 1 macro. 
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Macro Library name Macros defined 
“S$255SDUAZB:(SYSLIBISTARLET.MLB:2 0 (it(‘(Csé~;~S;*~™S 0 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRE SSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSPOWCC/OBJ=OBJ$:OTSPOWCC MSRC$:MTHJACKET/UPDATE=(ENHS :MTHJACKET) +MSRC 
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